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Abstract. In this article we derive some explicit decay estimates for the dual 
system of a basis of functions polynomially localized in space. 

1. Introduction 

A spline- type space 5 is a closed subspace of L'^{M.'^) possessing a Riesz basis 
of functions well localized is space. That is, there exists a family of functions 
^ S and constants < ^4 < i? < +oo such that 



(1) A\\c\\%<\\Y,c,M\l.<B\\c\\ 



holds for every c G and the functions {fk}f^ satisfy an spatial localization con- 
ditioiQ. 

In a spline-type space any function in f £ S has a unique expansion / = Ckfk- 
Moreover, the coefficients are given by Ck = ifigk), where {gk}k C S* is the dual 
basis, a set of functions characterized by the biorthogonality relation {g^, fj) = Sk.j- 

The general theory of localized frames (see [6], l5l and T) asserts that the 
functions forming the dual basis satisfy a similar spatial localization. This can 
be used to extend the expansion in ([T]) to other spaces, so that the family {/fc}^ 
becomes a Banach frame for an associated family of Banach spaces (see [4 and 
[6]). In the case of a spline- type space S, this means that the decay of a function in 
S can be characterized by the decay of its coefficients and that, in particular, the 
functions {fk}k form a so called p- Riesz basis for its L^-closed linear span, for the 
whole range 1 < p < oo. 

The purpose of this article is to derive, in some concrete case, explicit bounds 
for the localization of the dual basis. We will work with a set of functions satisfying 
a polynomial decay condition around a set of nodes forming a lattice. By a change 
of variables, we can assume that the lattice is 'E^. So, we will consider a set of 
functions {fk}^ ^ L'^{M.'^) satisfying the condition, 

(2) \fkix)\ <C {l + \x ~ kiy " , X eW'- andkeZ'^, 

for some constant C. This type of spatial localization is specifically covered by the 
results in [6], but the constants given there are not explicit. We will derive a poly- 
nomial decay condition for the dual basis {gk}^' giving explicit information on the 
resulting constants. This yields some qualitative information, like the dependence 
of theses constants on A, C and s and the corresponding p-Riesz basis bounds for 
the original basis. In particular, it implies that if a family of Riesz basic sequences 



^In particular, the series in the equation is required to converge unconditionally. 
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is given, where all the functions satisfy the concentration condition in equation 
([2]) uniformly, and have a uniform lower basis bound, then the corresponding dual 
systems are also uniformly concentrated in space. 

Since the localization condition in equation ([2]) is stable under small perturba- 
tions, the explicit estimates in Theorem ([T]) can be used to derive various kinds of 
stability conclusions. In particular, it could be used to derive jitter-error estimates 
for sampling in general spline-type spaces. 

The results in [6] prescribe polynomial decay estimates for the dual basis similar 
to those possessed by the original basis. As a trade-off for the explicit constants 
we will not obtain the full preservation of these decay conditions (see Remark |4]) . 
Nevertheless, any degree of polynomial decay on the dual system can be granted, 
provided that the original basis has sufficiently good decay. 

Finally observe that, although the basis {fk}^ is assumed to be concentrated 
around a lattice of nodes, the functions fk are not assumed to be shifts of a single 
function. In particular. Theorem [T] below allows for a basis of functions whose 
'optimal' concentration nodes do not form a lattice but are comparable to one - 
for example, it is possible to consider a system of functions constructed as integer 
shifts of a single one and then translate arbitrarily finitely many of them. The 
'eccentricity' of the configuration of concentration nodes is, however, penalized by 
the constants modelling the decay. 

2. Assumptions and statements 

Theorem 1. Let C > 1, and let t > d be integers. Let s > d + 1 be a real number. 
For k ^ li'^ , let fk : K"^ C be a measurable function such that 

\fk{x)\<C{l + \x-k\y\ {xeR''). 

Suppose that {fk}k ^■^ Riesz basis for S, its closed linear span within L^ , with 
bounds Q < A < B < oo (that is, equation ([T]) holds.) Then, the dual functions 
satisfy, 

\gk{x)\ < D {I + \x - k\y\ (xeR''). 

where D is given by, 

for some constant E > that only depends on the dimension d. 
Remark 1. The constant E can be explicitly determined from the proof. 

Of course, by a change of variables. Theorem [T] can be extended to a general 
lattice of the form PZ'^, where P in an invertible matrix. 

Since in the statement of the theorem we are allowing a constant that may 
depend on d, we can use any norm |-|. For convenience, we will use the max- norm 
\x\ = max{|a;i| , . . . , For the rest of the note, the functions fk^gk and the 

constants A, B,C, s,t will be fixed. 

3. Definitions and notation 

For a function /, the norm ||/|| (with no subscript), will denote its L^ norm. 
Given a matrix L = {lk,j)k,i£i'^ and 1 < /i < d, let 

Dh{L)k,j := {kh ~ jh)lk,j- 
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We denote by the norm of L as an operator on ^^(2"*). 

For two functions / and we write / < 5 if the exists a constant c > depending 
only on the dimension d such that / < eg. Furthermore, for functions /, g depending 
imphcitly on u > we say what / <„ g if there is a constant c > depending only 
on the dimension d such that / < c^g. 

Let us define the matrix M = {'rnkj)k.jeZ'^ by 

Wfej := ifkjj) ■ 

Since {fk}k is a Riesz sequence, AI is invertible. Moreover, ||M^^|| < A^^ and 
We also note the following estimates. The proof uses some techniques taken from 

Lemma 1. 

(a) Wu := ^ (1 + ^ ( 1 + ) ' keZ'^,u>d. 

(b) ^(l + |A:-j|)-«(l + |j|)-«<(l + |A;|)-", keZ'',u>d+l. 

j 

(c) / (1 + la; -2;|)-"(l + |y|)-"d2/< (1 + 1x1)-", x&R'',u>d+l. 



Remark 2. For items (h) and (e), the fact that u>d+\ is used to grant that the 
eonstant does not depend on u. 

Proof. First observe that for u > d, J^a + dx < ^j^- 

To prove (a) we simply estimate, 

E(i + i^ir"< E / ii + \k\r-dx. 



fcG 



For X G [0, 1]'^ + /c, 1 + |x| < 2 + l/cl < 2(1 + |fc|), so (1 + < 2"(1 + |a;|)- 

Therefore, 

E(i + ifcr"< / {i + \x\)-dx 

.alA k&-^ ifO.lJ'i + fe 

2" / (l+|a;|)-"da;<2"^— . 
Jgd u-d 



< 



So, for u < d + 1, 



u — d \ u — d 

k ^ 

For u> d+l, 

E(i + i^r" < E(i + « 1 < (1 + ^ 



u — d 
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For (b), we split the sum, 

^(l + |fc_j-|)-(<^+i)(l + |j|)-(rf+i) 
j 

b1>|fc|/2 

b1<|fe|/2 

For the first sum, if \j\ > \k\ /2, tlien 1 + |i| > 1 + \k\ /2 > 1/2(1 + and 
(1 + IjD-^'^+i) < 2(<*+i)(l + \k\)-^'^+'^\ For tlie second sum, if \j\ < \k\ /2, then 
\k-j\ > |A:|/2, and as before, {1 + \k - < 2(<'+i)(l + |A;|)-(''+i). 

Therefore, 

j 

\lil>|fc|/2 b-|<|fe|/2 
< (1 + \k\)-^''+^\ 

We now observe that, for any a > and all k 

{i + \k\r <{i + \k-j\rii + \j\r. 

Consequently, with a ^ u — d — 1, 

E(i +\k- + bi)-" (1 + < + 1^ - j\r^'^'\^ + b'i)"^'+'^ 

Therefore, 

E(i + |fc-il)-"(i + lil)-"<(i + |fc|)-"- 

This is the desired estimate. Assertion (c) follows similarly. □ 

4. Proofs 

We will first prove some lemmas and claims. 
Lemma 2. Let {cfc}^. C C 6e a sequence such that 

\ck\<a(l + \k-j\y\ keZ", 
holds for some constant a > and j ^IJ^ . Then the series defining the function, 

/ = E ^^^^ 
k 

is pointwise convergent and f satisfies 

\f{x)\<taC{l + \x-j\)-\ x&W' 



EXPLICIT LOCALIZATION ESTIMATES FOR SPLINE- TYPE SPACES 



5 



Proof. Given x e R'^, let ko e Z'^ be such that x e [0,1]'^ + kg. Then, for any 
k e i(l + |a; - < 1 + |/co - ^1 < 2(1 + |a; - A:|). Therefore, 

2-*(l + \x~ k\)-' < (1 + \ko - k\)-' < 2*(1 + |a; - 

Now we estimate, 

E ic^i 1^ < "^E(i + 1^ - (1 + - ^1)"' 

k k 

<aC^(l + |fc-jr*(l + |x-fc|)-* 
fe 

< aC2' ^(1 + |fc - (1 + |fco - fc|)~* 



<t aC{l + \ko - <t aC{l + \x- jI)"*, 

where we used that s > t and that t — d > 1. 

This proves that the series defining / is absolutely convergent at x and that / 
satisfies the required inequality. □ 

Claim 1. Let < u < t be an integer and 1 < h < d. Then, 

\\Dh{M)\\<C^Ws^t. 

Proof. For j, k ^1/^, first observe that, since s — d>t>l 

\{fk,f,)\<C^ I {l + \x-k\)-\l + \x-3\)-^dx<C\l + \k-j\)-\ 

Consequently, 

\DUM)k,,\<{i + \k~j\r\{fkj,)\ 

< c^i + \k~ j\Y{i + \k- = c^{i + \k- 

Therefore 

supJ2\Dl{M)k,,\ , sup^ \DUM)k,,\ < C^Ws-t. 
3 ' k 

Consequently, by Schur's Lemma 

\\Dl{M)\\<C^W,^t. 

□ 

Claim 2. Let < u < t be an integer and 1 < h < d. Then, Z))^(M^^) is a 
bounded operator on 

Proof. In Claim [T] we have shown that the matrix M satisfies the off-diagonal decay 
condition, 

< A'(l + lfc-J|)-^ 

for some constant K > Q that depends on d and C . 

Jaffard's Theorem ([5], see also [B]) implies that satisfies a similar off- 

diagonal decay condition, 

M-] <x'(l + |fc-iI)-^ 
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for some constant K' > 0. A calculation similar to that in Claim [U shows that for 
< u < t, Z3Jj'(A/^^) is a bounded operator (and gives a norm estimate depending 
on the unknown constant K'). □ 

Claim 3. Let I < h < d. Then, 

Proof. Let u, 1 < u < t be an integer. Consider the map L i— > Dh{L), defined for 
those matrices L such that Dh{L) is a bounded operator. We know from Claims [T] 
and m that we can apply Dh to both M and A/^^, up to t times. The map Dh is 
a derivation from its domain into B{£^). This means that it satisfies the equation 
Dh{AB) = Dh{A)B + ADh{B), provided that DhiA) and DhiB) are both weU 
defined. Iterating this equation yields, 

(3) = Dl{M-'M) = J2 (^^D\XM-')Dl-\M). 

Since all the operators involved in the last formula are bounded, we can associate 
factors to obtain, 

(4) Dh{M-^) = (''^Di{M-^)Dl-\M)M-\ 

1=0 ^ ^ 

It follows that 

U-l , s 

1=0 ^ ^ 

Consider the numbers w„ := maxo<„'<u||D^ (Af^^)||. Using Claim [TJ we ob- 
serve that A < inf {|A| : A G cr(M)} < ||M|| < C'^Ws-t, where a{M) denotes the 
spectrum of the Gramian matrix AI as an operator on 

Consequently 1 < A'^C^Ws-t and the numbers Vu satisfy, 

U-l , s 

(=0 ^ ' 

va < A-\ 
Iterating the last relation yields. 



□ 



Now we can prove the main result. 



Proof of Theorem\^ For G if", let Ckj :— {gk,gj). Since M/^j^ = Ckj, using 
Claim |3] we get, 

M < \\Di{M-')e,h <t^ A-('+'^C''Wl_,. 
Recalling that we are using the max-norm, we obtain 

For k 7^ j, \k — j\ > 1, so (1 + |/c ~ <t \k — jf. Consequently, 



' s-f 



EXPLICIT LOCALIZATION ESTIMATES FOR SPLINE- TYPE SPACES 7 

holds for k / j. For k = j we use that |cfcfc| — WgkW^ < A^^ . Indeed, since {/fc}^ 
is a Riesz sequence with lower bound A, it is also a frame sequence with the same 
lower bound. That is, the inequality, J^k \{fJk)\^ > A\\ff holds for all f G S. 
Setting / — gj, it follows that WgjW^ < A~^. 
Since Wg , C > 1 and we have that 

(1 + \k-j\y \ck.j\ <t^ + A-')C^'Wl,. 

for aU fc, j e Z'^. 

Since for each k E 1,'^, gi^ = Ck.jfj: Lemma [3] implies that 

\gu{x)\ <t2 + + |a; - x e M'^. 

By Lemmam W^^t < (l + (^rbi) ) ■ Therefore Wl_t <t (l + (jriZTiy ) *■ □ 
We conclude that, 

(5) 

|5fc(x)| <t2 + {l + is-t- (1 + |x - fc|)-*, a; £ M'^. 

This is almost the desired estimate; there is an extra A^^ term. We get rid of it 
using the following argument. 

Let a be a positive real number. Consider the functions fk '■— afk- They 
form a Riesz sequence with lower and upper bounds a^A and a^B respectively. 
Moreover, the functions forming the dual basis are given by gk a~^gk. The 
family | /fe | satisfies a concentration condition similar to the one possessed by the 

original functions, but with a constant aC. We apply the estimate in equation ^ 
to this new family and learn that, for a; G M'*, 

\9k{x)\ <t2 (a(-2t-2)^-(t+i) ^ ^-2^-1^ ^(2t+i)^2t+i (1 + _ ^ _ ^yiy (1 + 1^ _ k\)-'. 

Therefore, for a; G K'', 

|.gfc(:E)| <t2 + (i + (5 _ t _ {I + \x - fc|)-*. 

Letting a — > 0+, we get the desired estimate. 

5. Final remarks 

Remark 3. The most delicate part of the proof was the justification of the formal 
computations in Claim\^ that allowed us to solve \\D'^{M~^)\\ recursively from 
the binomial formula. In order to associate factors, we needed to know that 
belongs to the domain of D"^ . 

To see why this is important, let us consider the case where the original basis is 
formed by the integer shifts of a single function. In this case, the Gramian matrix 
M is a convolution operator, having some sequence a as kernel. The matrix 
is also a convolution operator and has a kernel b that satisfies, 

(6) a* b = S. 

The decay of a and b can be reformulated in terms of smoothness estimates for their 
Fourier transforms a and b. The argument in the proof of Theorem{I\ amounts, in 
this case, to transferring smoothness estimates from a to its pointwise inverse b by 
an iterated application of the Leibniz product rule ( cf. Equation ^ .) 
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The obstacle to derive Equation formally from Equation O is that the latter 
equation does not determine, by itself, the derivatives of M^^. For example, when 
a is a finitely supported sequence, Equation ^ is a recurrence equation in b, that 
has many solutions even if a has no zeros. The sequence b that we are looking for 
(that is, the kernel of M^^) can be singled out as the only solution of Equation 1^ 
that belongs to . 

In the case of a basis formed by the integer shifts of a single function, the jus- 
tification we need follows from some careful regularization argument for Sobolev 
spaces. In our case, this justification was done in Claim [H by resorting to Jaf- 
fard's result Another possible approach would be to use the general theory of 
unbounded derivations, in particular the results in [2] and [9] . However, this would 
require adapting those results to non-densely defined derivations. 

The use of derivations is somehow implicit in Jaffard's original proof 8 . A re- 
cent article from Grochenig and Klotz [J exploits in great depth the use of deriva- 
tions to study off-diagonal decay of infinite matrices and approximation by classes 
of banded matrices. 

Remark 4. As observed in the introduction, in Theorem [H the decay condition 
on the original basis {fk\k ^'^^ shown to be fully shared by the dual basis {gk}k 
(although the results in [6J show that the full decay condition is actually preserved.) 
This is due to the kind of objects used to bound the decay of the entries of M and 
M^^ . According to Remark\^ above, in the case of a basis formed by the integer 
translates of a single function, the estimates given amount to smoothness estimates 
for the symbol t of some convolution operator. In Claims\^ and\^ we bounded the 
size of the entries of a matrix by means of its P' operator norm and controlled 

that norm by interpolating its £^ — s- £^ and £°° i°° norms (by Schur's lemma.) 
This would correspond in the case of a convolution operator to bounding the L°° 
norm of its symbol t, from above by its norm and from below by its norm. 
This accounts, in that case, for the loss of some precision in the estimates. 

Remark 5. Finally, it should be noticed that the techniques used in this article 
could be applied to a general polynomially self-localized basis in the abstract sense 
of [5], where a basis is called self-localized if its Gramian matrix presents certain 
off- diagonal decay. 
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